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Abstract 

A summary of discussions on selected topics related to QCD phase diagram, phase transition, 
critical point, fluctuations and correlations at the Quark Matter 2009 conference are presented. 



1. Introduction 

Understanding the properties of the nuclear matter when subjected to extremes of temperature 
and density is one of the goals of the relativistic heavy-ion collision program. The knowledge of 
the properties manifested in these collisions should have some connection to our understanding 
of the early universe and how hadrons acquire their masses. QCD suggests we should at least 
expect two types of transitions for nuclear matter in limits of high temperature (T) and densities : 
deconfinement transition and chiral phase transition. Theoretically the deconfinement measure is 
the order parameter, Polyakov loop (L(T) ~ lim r ^ M exp{- V(r)/T), where V(r) is the potential 
between a static quark-antiquarkpair separated by a distance r). The measure of chiral transition 
is the order parameter, the chiral condensate (tfrt//)(T). For a confined phase for T < T c (critical 
temperature for transition) we have L(T) and (tfrt//)(T) + 0, and for a deconfined phase for 
T > T c with L(T) + and (tfnfr)(T) 0. In this proceedings we summarize the theoretical and 
experimental understanding of the QCD phase diagram as was discussed at the QM2009. 

2. Lattice QCD results 

Lattice calculations simulate a quantum statistical ensemble in thermal equilibrium at fixed 
T with partition function Z = Tr e ~ H / T t where H is the QCD Hamiltonian. The lattice spacing 
a and N T sites in imaginary time are related to T as aN, = l/T. Usually a is varied to change T 
keeping N T fixed. An approach based on "T integral method" was presented, where N T is varied 
and a fixed to change T QJJ] . There are few aspects to be noted while interpreting lattice results. 
To connect with reality, lattice results needs to be presented as extrapolations to continuum limit. 
The calculations are now done on lines of constant physics, along which the variations in ob- 
servables can be attributed to changes in temperature and not also to changes in the Hamiltonian. 
Choice of proper actions (see recent work |2[]), spatial volume and setting of quark masses are 
important for interpreting the results. Results are also dependent on the choice of number of 
quark flavours. The current view is, computation with 2 light quarks and a heavy strange quark 
(2+1) are close to a realistic picture. At the conference, several studies from Lattice QCD were 
presented. These provided insights on the order of phase transition, transition temperature (T c ) , 
Equation-Of-State (EOS) and QCD critical point (QCP), some of these are discussed below. 
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Figure 1 : Left: Susceptibility for the light quarks for 
N T =6 as a function of 6/g 2 , where g is the gauge cou- 
~ pling (T grows with 6/g 2 ) (J. Right: Derivative of 
logarithm of partition function in canonical ensem- 
-j 2 ble formalism as a function of the quark number den- 



sity d. 



2.1. Order of phase transition 

QCD calculations on lattice at high temperature and ps — MeV has established the quark- 
hadron transition to be a cross-over yfl. FigureQ]shows the lattice chiral susceptibility x(N s ,N t ) 
= d 2 /(dm 2 d )(T /V)- logZ, where m ut / is the mass of the light u,d quarks, N s is the spatial extension, 
N T euclidean time extension, and V the system volume. The susceptibility plotted as a function 
of 6/g 2 (g is the gauge coupling and T grows with 6/g 2 ) shows a pronounced peak around the 
T c . The peak and width are independent of volume (varied by a factor 8) thereby establishing 
the transition to be an analytic cross-over |3j]- For a first-order phase transition the height of the 
susceptibility peak should have been oc V and the width of the peak oc 1/y, while for a second- 
order transition the singular behaviour should have been oc V a , a is a critical exponent. A lattice 
result discussing the nature of phase transition at finite T and baryon density J4J] was presented. 
It was based on idea that the chemical potential (p*) which gives a minimum of the effective 
potential does not increase when it passes through the mixed phase. Whereas for a cross-over 
transition it shows a monotonic increase with density when T > T c . The results are shown in 
the right panel of Fig. Q] This calculation which uses P4-improved staggered quark action, in 
canonical ensemble with m n = 700 MeV and number of quark flavour of 2, results in a first order 
phase transition for T/T c < 0.83 and (i q /T > 2.3 0]. 

2.2. Transition Temperature 

The point of sharpest change in temperature dependence of the chiral susceptibility (x^), the 
strange quark number susceptibility (x s ) and the renormalized Polyakov-loop (L) are used to es- 
timate the QCD transition temperature in lattice calculations. The recent results on chiral and de- 
confinement phase transition temperatures were presented by two groups (HotQCD/RBC [5] and 
Budapest/Wuppertal [6]). Using the observable;^ /T 2 the HotQCD/RBC and Budapest/Wuppertal 
groups get the chiral phase transition temperatures as 192 (4)(7) MeV and 152(3)(3) MeV respec- 
tively. Using the observable L the deconfinement transition temperature remained the same for 
HotQCD/RBC and is 170(4)(3) MeV from Budapest/Wuppertal group. The two groups differ in 
the transition temperatue for both chiral and de-confinement transitions. The possible sources of 
differences in the two approaches could be ambiguity in locating T c for a cross-over, physical 
observable used to set the scale, preferred renormalization of chiral susceptibilities and choice of 
actions. It was suggested at the conference to study the QCD thermodynamics with a theoreti- 
cally firmly established Wilson type fermion discretization. Such a large difference in the value 
of T c has serious consequences for heavy-ion phenomenological studies. As energy density (e) 
~ T 4 , the lack of accurate determination of T c leads to about 60% uncertainty in e at T c . 
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Figure 2: Left: Results for energy den- 
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2.3. Equation of State 

Figure [2] (left) shows the variation of e/T 4 vs. T J3, Si- As e/T 4 oc effective degrees of 
freedom (dof) of the system (g e frX the trend shows a change in g e g below and above T c . This 
is interpreated as a transition from a state with hadronic dof to a state with quark-gluon dof. 
At high T the e/T 4 is about 10-15% lower compared to Stefan-Boltzman (SB) limit indicating 
that the matter still has strong interaction effects. The deviations from ideal gas behaviour could 
be understood in terms of effective thermal masses of quarks and gluons, and a consequence 
of this sizable interaction at high T could lead to the existence of coloured resonance states. It 
was shown that at very high T ~ 10 5 - 10 7 MeV the lattice calculations approach the SB limit 
and has good agreement with calculations from perturbation theory [9]. The functional relation 
between pressure (p) and e is the EOS. The velocity of sound in the medium (c 2 s ) is dp/de at 
fixed entropy, p/e is plotted as a function of e 1 ' 4 computed from Lattice in Fig. [2] (right) fl. 
This is an important ingredient to heavy-ion phenomenology calculations, as the c s decides the 
rate of cooling of the high-r system and provides input on the system composition at different 
evolution times. For the dependence of thermodynamics of QCD plasma on number of colors 
see Jioll . 



3. Experimental results on Charge Correlations 

The STAR experiment at RHIC presented preliminary results on charged hadron azimuthal 
correlations based on 3-particle correlation technique II in . The results from Au+Au collisions 
at yfsim = 200 GeV at midrapidity for 0.15 < pj < 2 GeV/c between same charged and oppo- 
site charged hadrons with respect to reaction plane (plays the role of third particle) are shown in 
Fig. |3] The observable, (cos(0 fl + <pp- 2^^)) represented the difference between azimuthal cor- 
relations projected onto the direction of the angular momentum vector and correlations projected 
onto the collision event plane. The difference between the same charge and opposite charge cor- 
relations could not be explained by models such as HIJING and UrQMD and by incorporating 
realistic values for the elliptic flow in such simulations. The signal seems to be consistent with the 
predictions for existence of metastable domains in QCD vacuum leading to local Parity violation 



in the vicinity of deconfinement transition expected to be achieved in heavy-ion collisions 11211 . 
For such a phenomena where the massless quarks can change their chirality due to interactions 
with gluon fields, there could be separation of positive charges from negative charges along the 
direction of angular momentum of the collision as a result of large magnetic fields reached in 
the collisions (especially in non-central collisions). Deconfinement allows for the possibility of 
quarks traveling over distances greater than nucleonic scales and chiral symmetry restoration is 
essential, because a chiral condensate will tend to erase any asymmetry between the number of 
right- and left-handed fermions. The observable presented is parity-even 111 ill , making it suc- 
ceptible to physical processes not related to parity violation effects. Thus understanding of the 
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Figure 3: Azimuthal charge correlations in data com- 
pared to simulation results for 200 GeV Au+Au flUl . 
Blue symbols mark opposite-charge correlations, and 
red are same-charge. Markers connected by solid lines 
represent the data. 



physical background is crucial. The signal was also observed to be not restricted in the low pj 
region as naively might be expected for parity-violation effects. 

4. QCD Critical Point 

4.1. Theory Calculations 

Several QCD based models predict the existence of an end point at high ps for the first order 
phase transition in the QCD phase diagram [1X311 . However the exact location depends on the 
model assumptions used. A summary of these model results can be found in Ref. ffl . One 
such calculation based on linear <x model coupled to two flavors of identical mass quarks was 



presented at the conference [15]. This work extends the previous investigations by including 
thermal fluctuations of the meson and fermion fields. It is observed that when all other parame- 
ters of the model are fixed, the calculation allows for existence of zero, one or two critical points 
in phase diagram depending on the value of the vacuum pion mass. A typical result from this cal- 
culation for varying mass of pion is shown in Fig. [4] Given the ambiguity in predictions of QCP 
in models, studies on lattice was expected to provide reliable estimates. However lattice calcula- 
tions at finite ps have important issues to be addressed. Typically, for any lattice computation one 

needs to evaluate the expectation value of an observable X, (X(m,,)) = ^°^* P ^ S , G ^"'^ Dt * M< -'"'' > , 

r J DUexp(-S G ) Det M(m s ) ' 

where M is the Dirac matrix in x, colour, spin, flavour space for sea quarks of mass m s , So is 
the gluonic action, and the observable X may contain fermion propagators of mass m v . The 
Det M for non zero p is not positive definite, hence numerical methods of evaluation of the ex- 
pectation values is difficult, this is commonly refered to as the sign problem. There are several 
ways suggested to overcome this issue, (i) Reweighting the partition function in the vicinity of 
transition temperature and p — [16], (ii) Taylor expansion of thermodynamic observables in 
p/T about p = B17I1 and (iii) Choosing the chemical potential to be imaginary will make the 



ferminonic determinant positive II 1 811 - The first two methodologies yield an existence of QCP 



(as shown in Fig. |4|, whereas the thrid procedure gives a QCP only when the first co-efficient 
in the Taylor expansion of generic quark mass on the chiral critical surface (m c ) as a function of 

f 1 ^ (/nf}) = 1 + Tjk=\ c k {yr) -* * s positive. The calculations seem to support a negative value 
of c\ indicating an absence of qcp3 This result was intensely discussed at the conference with 



Developments after QM2009 have shown that with decreasing lattice spacing, c\ becomes positive and the expected 
picture of the phase diagram seems to be recovered (O. Philipsen at CPOD2009). 
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Figure 4: QCD phase diagram in T vs. /jg plane showing the QCP from various model calculations. For details see text. 
From left to right the plots are adapted from Refs. fT3l . ll3l . and fl^l . respectively. 

suggestion to do the calculation with larger spatial volume, check the stability of the results for 
different values of N T and for a realistic 2+1 flavour in continum limit. A new calculation on 
lattice using the Canonical ensemble was presented at the conference showing the existence of 
an end point for Te ~ 160 MeV and fj.£ ~ 600 MeV lEoll . Further improvements are expected 
when calculations will use a realistic pion mass, instead of m n = 700 MeV used in the presented 
work. 

4.2. Signatures of QCD Critical Point and Experimental Results 

The characteristic signature of QCP is large fluctuations in event-by-event conserved quan- 
tities like net-charge, net-baryon number and net-strangeness Izill . The variance of these dis- 
tributions ({{5N) 2 )) are proportional to square of the correlation length (£). The finite size and 
finite time effects attained in high energy heavy-ion collisions, limits the value of the £ achieved 
in the collisions. Model calculations suggest it could be small (2-3 fm) 1E2I1 . thereby making it 
extremely challenging to measure in the experiments. Motivated by the fact that non-Gaussian 
features in above observables increase if the system freezes-out closer to QCP, it has been sug- 
gested to measure higher moments (non-zero skewness and kurtosis indicates non-Gaussianity) 
of net-charge or net-baryon number distributions. Further it has been shown that higher mo- 
ments (((5N) 3 ) ~ 5 and ((6N) 4 ) ~ £ 7 ) have stronger dependence on £ compared to variance and 
hence have higher sensitivity 12311 . Another important reason to look for higher moments comes 
from lattice calculations of the quadratic and quartic net-charge, net-baryon and net-strangeness 
(these denoted by Nx) fluctuations (evaluated at vanishing chemical potential), x\ — Wf(^x) 
and X4 = yfi ((N x ) - 3(A^) 2 ) 12411 . x\ show a rapid rise in the transition region whereas fhe^ 
show a maximum at T c . This maximum is most pronounced for the baryon number fluctuations 
and would diverge at QCP due to long-range correlations |19]. So studying higher moments 
allows for a connection between QCD calculations on lattice and experimental data. Experimen- 
tally it is difficult to measure all the produced baryons. This issue was addressed theoretically, 
where it was shown that net-proton number fluctuations would faithfully reflect the singularity of 
the charge and baryon number susceptibility due to the iso-spin blindness of a field [23[]. Hence 
looking for non-monotonic variation of higher moments of net-proton distribution as a function 
of yfsim (or T - fig) is sufficient to locate the QCP. The first results on moments of net-proton 
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Figure 5: Moments of net-protons vs. centrality (left panel) 12511 . slope of anti-baryon to baryon ratio vs. pj (middle 
panel), system size dependence of mean pj fluctuations (right panel) 1 26]. See text for details. 



distribution at RHIC energies was presented at the conference and shown in Fig.[3][25]. The evo- 
lution of the moments from preipheral to central collisions follow the expectations (dashed lines) 
from central limit theorem (equations shown in Fig.|5]l. This study at fig < 30 MeV provides an 
understanding and a formulation for the physics background for the observable expected to be 
sensitive to long range fluctuations as expected for QCP. 

Several other interesting signatures of QCP and associated results from SPS were discussed 
at the conference, one of them was based on the idea that the presence of a critical point deforms 
the isoentropic trajectories in the T - ps phase plane. The critical point serves as an attractor 
of the hydrodynamical trajectories. This later feature was debated at the conference. If such a 
scenario exists, it would lead to an experimental signal of drop in pi p (~ e -2 ^/ 7 ") vs. pj (in 
intermediate region) i27ll . The suggestion was put to test with the existing data, by looking at the 
slope of pip vs. pj at available yfsim at SPS and RHIC energies. As shown in the Fig.[5](mid- 
dle panel) there is no large drop in slope is observed for intermediate pj range. Another novel 
idea discussed was based on study of critical dynamics around QCP with relativisic dissipative 
hydrodynamics M28I1 . Using the idea of coupling the density fluctuations to thermal energy it was 
shown that sound modes around QCP will be suppressed. This will then lead to disappearance 
of mach-cone like signals observed at top RHIC energies J29I1 . NA49 experiment presented mul- 
tiplicity and mean pj fluctuation results as a function of beam energy and ion size B26I1 . Within 
the experimental acceptance no strong non-monotonic dependence of fluctuations was observed 
for Pb+Pb collisions as a function of yfsim (or pg)- However the system size dependence (using 
the central C+C, Si+Si and Pb+Pb data) seems to indicate a non-monotonic behaviour around 
Pb ~ 250 MeV and T ~ 178 MeV. This was shown to be consistent with a theory calculation 
including a QCP and with £ ~ 6 fm and folding in the experimental acceptance (shown as solid 
line in Fig.[5]nght)). This study further underscores the need for a more detailed investigation in 
the upcoming critical point search programs at RHIC and SPS. 



5. Phase Diagram 

The current understanding of the QCD phase digram as discussed at the conference is de- 
picted in Fig. [6] From the QCD calculations on lattice it is now established theoretically that the 
quark-hadron transition at ps — MeV is a cross-over 0], The critical temperature for a quark- 
hadron phase transition lies within a range of 170-190 MeV (vertical band in Fig. H H. 
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Figure 6: Current understanding of T-fig phase diagram for nulear matter from theoretical QCD calculations, experi- 
mental results and discussions at the conference. RHIC QCP search: ysjviV ~ 39 - 5 GeV (collider); SHINE: lighter 
nuclei, ^snn ~ 17-5 GeV (fixed target) and CBM: high luminosity, yiwiv ~ 10 - 4 GeV (fixed target). 



Most calculations on lattice also indicate the existence of QCD critical point for pg > 160 
MeV EEllMl. The exact location is not yet known unambi giou sly. Two such predictions 
computed on lattice are shown in Fig. for a T c of 176 MeV [M d- High energy heavy- 
ion collision experiments have seen distinct signatures which suggest that the relevant degrees 
of freedom at top RHIC [30] and SPS energies IB ill in the initial stages of the collisions are 
quark and gluons. Specifically the direct photon data is used to show that the initial temperature 
{T initial, Fig.|6]l reached at RHIC and SPS is greater than T c predicted by lattice 1 32] . Further the 
understanding of suppression in high pj hadron production in heavy-ion collisions relative to 
p+p collisions at RHIC requires a medium energy density > 1 GeV/fm 3 (critical energy density 
from lattice for a phase transition) H30TI . The experiments have also measured the temperature at 
which the inelastic collisions ceases (Chemical freeze-out) and ellastic collisions ceases (Kinetic 
freeze-out) [33]. These temperatures (as shown in Fig. [6j are extracted from the measured par- 
ticle ratios and transverse momentum distributions using model calculations which assume the 
system is in chemical and thermal equilibrium. It is interesting to note that the difference between 
the two freeze-out temperatures becomes smaller at high pg (estimated at chemical freeze-out). 
New experimental programs at RHIC Q, SPS JH], FAIR JH] and NIC A facilities have been 
designed to search for the QCD critical point in coming years. Whereas the experimental pro- 
gram at LHC (probing the region of pg ~ MeV of the phase diagram) will provide an unique 
opportunity to understand the properties of matter governed by quark-gluon degrees of freedom 
at unprecedented high initial temperatures (higher plasma life time) achieved in the Pb+Pb colli- 
sions at 5.5 TeV [37]. A novel theoretical proposal was made at the conference on the existence 
of a quarkyonic phase around pg values corresponding to AGS energies [38]. This is in addi- 
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tion to confined and de-confined phases. The matter in such a phase is expected to have energy 
density and pressure that of a gas of quarks, and yet be confined. Baryon-Baryon correlations to 
look for nucleation of baryon rich bubbles surrounded by baryon free regions was discussed as a 
signature of such a phase i39ll . 
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